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We formulate the statistical mechanical description of liquid systems for both polarizable and 
polar systems in an electric field in the E-ensemble, which is the pendant to the thermodynamic 
description in terms of the free energy at constant potential. The contribution of the electric field 
to the configurational integral Qn{E) in the E-ensemble is given in an exact form as a factor in 
the integrand of (5jv(E). We calculate the contribution of the electric field to the Ornstein-Zernike 
formula for the scattering function in the E-ensemble. As an application we determine the field 
induced shift of the critical temperature for polarizable and polar liquids, and show that the shift 
is upward for polarizable liquids and downward for polar liquids. 
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I. INTRODUCTION 



The behavior of liquid systems in an external electric 
field has been the subject of experimental and theoreti- 
cal interest over many decades [H-[20|- ^ recent example 
are block copolymer melts which are of particular inter- 
est due to field induced uniform macroscopic alignment 
in the microphasc separated state (and references 

therein), which is of basic importance for applications 
using self assembled block copolymer structures for pat- 
terning and templating of nanostructures [2l| . A more 
basic question concerns the shift of the critical temper- 
ature for low molecular weight systems in an electric 
field. Although this issue has been the subject of in- 
vestigations over many years ^ microscopic ex- 
planation of the shift has not been provided. Most of 
the existing theoretical work on the contribution of the 
electric field to thermodynamic quantities is based on a 
phenomenological description in the framework of macro- 
scopic electrodynamics. It is well-known 4] that for ther- 
modynamics in an electric field one should distinguish 
between thermodynamic potentials at constant dielectric 
displacement or constant charges F{T, V, D) and at con- 
stant electric field or potential F{T,V,'E). However, the 
available theoretical studies do not provide the tools for 
calculating the thermodynamic quantities and the corre- 
lation functions (in particular the scattering function) in 
the ensemble at constant electric field, which is of pri- 
mary experimental relevance. In this Letter we address 
this problem and establish the statistical mechanical de- 
scription of liquid systems in electric fields at constant 
electric field F{T, V, E) by introducing the T - F - E- 
ensemble following the analogy between the T — V and 
T—p ensembles of Statistical Mechanics [l^. The contri- 
bution of the electric field to the configurational integral 
Qn(E) in the T — V — E-ensemble is given in an ex- 
act form as a factor in the integrand of Qn^E). This 
E-cnscmble provides a systematic way for the calcula- 
tion of thermodynamic quantities, correlation functions 
in electric fields, and among others the calculation of the 
dielectric constant. We also derive the contribution of 



the electric field to the Ornstein-Zernike expression of 
the scattering function to quadratic order in atomic po- 
larizabilities in both Eg- and E-ensemble. Using these 
results we consider the shift of the critical temperature 
for both polarizable and polar liquids in an electric field, 
and give an explanation for the different sign of the shift 
in these systems. 

The article is organized as follows. Section |TT3 intro- 
duces to the formalism of the description of simple liquids 
in an electric field. Section IIIBI generalizes the formal- 
ism to polar systems. Section IIIII considers the shift of 
the critical temperature in simple and polar liquids, and 
the comparison with experimental results. 



II. STATISTICAL MECHANICS OF SIMPLE 
LIQUIDS IN AN ELECTRIC FIELD 

A. Description of polarizable systems in the 
E-ensemble 

Let us first review the statistical mechanics of a sys- 
tem of interacting molecules without permanent electric 
dipole moments with Hamiltonian Hq in an external elec- 
tric field Eq. The induced dipole moment of the ith 
molecule for a system of molecules in an external elec- 
tric field Eq is given by Q 



(1) 



where the summation occurs over all molecules, and 
the molecular polarization tensor a is defined by pf = 
afE^iji) with a"'3 = aj"l^. The tensor of dipole- 
dipole interactions T"^ is given by 



J 3 , 
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It is assumed in ([T]) and in the following that the diagonal 
elements Tu are zero. The total microscopic strength of 
the electric field at the position r is given by the sum 
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of the external field Eo(r) and the field of the induced 
dipoles Pj 



i?"(r)=i?«(r)-^r"'3(r-r,)p; 



(2) 



The total interaction energy of the induced dipoles is 
given by Q 

Hpoi = -i5]i?o(r-.)[l + «.T]->«,£;o(r,) 



(3) 



The last line in ^ corresponds to the electric energy of 
the induced dipole moments pf given by Eq. ([T]) in the 
external electric field i?Q (r^). The expansion of Hpoi in 
Eq. ^ to linear order in T yields 

Hpoi -'IY. ^oaEo + ^oaiT(r, - rj)ajEo + ■■■ . 

(4) 

The Cartesian indices are suppressed in The config- 
uration integral in the electric field is given by 

QNiEo) = / dVi • • • / d^r^e-^^"-^^-- , (5) 



where Hq is the Hamiltonian which contains the repulsive 
hard-core and the attractive van der Waals interactions 
between the molecules. Eq. ((S]) allows the statistical me- 
chanical computation of the free energy in a constant 
external electric field Eq which can be identified with the 
dielectric displacement D, because both obey the same 
Maxwell equation divEg — A-np^xt- 

The average of Hpoi given by Eq. @ over positions of 
molecules results in 
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(6) 



To calculate the integral over q in the first line 
of Eq. ([6]) we replaced the delta function (5'-'^''(q) 
by the normalized bell-shaped function (5a (q) ~ 
(f /27ra^)^/^ exp(-(l/2a^)q^). The free energy is related 
to Qjv by = Fo — ksTlnQN/V^, and does not contain 
the energy of the external electric field. The expansion of 
Qat in ([5]) up to linear terms in dipole-dipole interactions 
and the use of ^ results in the following expression of 
the free energy, which is exact to 2nd order in as 



Stt N 



(7) 



In thermodynamics the free energies at constant di- 
electric displacement F{T, V, D) and at constant electric 



field F{T, V, E) are related by a Legendre transform as 
F = F - EDF/47r with the differentials given by [l] 



dF = -SdT - pdV + — EdD, 
47r 

dP = ^SdT - pdV - — DdE, 
47r 



(8) 



where for simplicity homogeneous fields are assumed in 
the above expressions. 

We will now consider the question of the statistical me- 
chanical calculation of the thermodynamic quantities in 
the ensemble at constant potential. According to the gen- 
eral principles of construction of different ensembles in 
Statistical Mechanics [22] the relation between F and F 
should correspond to the transition from the T — V — Eq- 
ensemble to & T — V — E-ensemble. Following the known 
examples in Statistical Mechanics ■22||, as for example the 
relation between the T — V and T — p ensembles, we de- 
fine the T — V — E ensemble by the following expression 
for the configurational integral 



Qjv(E) 



L>Eo(r)exp [ (3 d^r 



EE(] 

47r 



X exp 



Qn{^q), 



(9) 



where the term ~ Eq corresponds to the energy of the 
external electric field. The integration over the field 
strength in Q occurs at every r, i.e. ([9|) is a functional 
integral. 

Indeed, the functional integral in ([9]) is Gaussian, and 
consequently the integration over Eq (r) can be performed 
exactly using the quadratic complement. We obtain 



QAr(E) = / drexp(-/3Fo)exp 



-ln(— det^-i) 



X exp 



^ I d\ I d^r'E{r)A-\r,r')E{r') 



(10) 



where dT denotes here integrations over positions 
ri, • • • , rjv of the particles, and the matrix A is defined 
by 



A{r, r') = S{r - r') - Aim{r){I + aTn)^ l,c 



(11) 



where / 5{r — r') is the identity matrix, and n{r) is the 
microscopic density n(r) = '^5{r — ri). The Cartesian 

i 

indices of A are suppressed. The expansion of A^^{r, r') 
under the integral in the last line of (fTU|) to the sec- 
ond order in atomic polarizability gives rise to a factor 
exp{-(3W2,poi) with 



= -lY.E'^af[AnS0'^d{r,-rj) 
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r,)] i?^. 



(12) 



The first term in p2p originates from interactions of in- 
duced dipoles with the electric field to the 2nd order. 
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while the second term is the 1st order contribution of 
the dipole-dipole interactions. For homogeneous field the 
term linear in as does not depend on the positions of the 
molecules, and consequently does not contribute to ob- 
servables. Eq. (fTO)) is the starting point for a calculation 
of the dielectric constant by summations of subseries of 
perturbation expansions in powers of the atomic polariz- 
abilities. It is easy to see that the preaveraging of A(r, r') 
in the exact expression (fTTj) according to A{r, r') 

{A{r, r')) = S{r - r') (l - Anna {I + 87ran/3)~^) results 

in the following expression for Qn(E) — exp{—(3F) = 
exp {-f3Fo + P J eE^(r)d^r/87r) with the dielectric con- 
stant 1/e = 1 — 4:T:na/ (1 + 87rna/3), which is equivalent 
to the Clausius-Mossotti relation. Effects of the electric 
field beyond can be also studied using (fTO|) . The Indet 
term in (fTU|) does not depend on the electric field but on 
the positions of the molecules. It has the following in- 
terpretation: Induced dipoles appear also in the absence 
of an external electric field due to thermal fluctuations. 
Their contribution to the interaction energy between the 
molecules can be absorbed into interactions between the 
molecules which arc independent of the external electric 
field, and are out of interest in the present work. 

We now will establish the contribution of the electric 
field to the static scattering function in the Eo-ensemble. 
We start with the Ornstein-Zernike equation ^] 



r2) = c(ri - r2) + n J (fr^ciri - r3)/i(r3 - r2). 



which expresses the correlation function h{ri — = 
g{ri — r2) — 1 through the direct correlation function 
c(ri — r2), where n = N/V is the average density, and 
the distribution function g{ri — r2) is defined by 



- r2) 



Jd^ri---Jd^rNe-l^"o-'^"p'" 



(13) 



In the case of homogeneous field we obtain from ([13]) to 
the lowest order in dipole-dipole interactions 

g{ri - r2) « go(i"i - ^2) - /3SoaT(ri - r2)aEo. 

The scattering function is expressed through the Fourier 
transform of the direct correlation function as 
5^^(q) = 1 — nc(q). Using the above relations and the 
Ornstein-Zernike formula we obtain the scattering func- 
tion in the electric field Eq in the vicinity of the critical 
point by taking into account the interactions of induced 
dipoles to lowest order as 



5-i(q)=r + cg2+47r/3n«2Ml) 



(14) 



where t — T — T^. 

To establish the contribution of the electric field to the 
scattering function in the E-ensemble we start with the 



expression of the radial distribution function 

V /■„_,, f„ f E(r)E„(r) 



47r 



d\^... / d\Ne'^"^'-^"-"\ (15) 



Qjv(E) 
Eg(r) 
Stt 



which is in accordance with the definition of the partition 
function in Eq. Carrying out the functional integra- 
tion over Eo(r) similarly to Eq. PU)) we arrive at 



5(ri,r2) 



1 



1/2 / 



where (5a'(E) and A{r^r') are defined by Eqs. (jlOp and 
(fTTj) . respectively. The lowest-order contribution to the 
direct correlation function from the electric field can be 
obtained by expanding A^^(r, r') in Eq. (jl6p similar to 
the corresponding expansion of Eq. (jlOp . which results in 



5(ri,r2) ~5o(ri - r2) exp (-^1^2,^0!) , (17) 

where the effective binary energy W2.P0I is given by 
Eq. (fT^ . Following the same steps as in the above deriva- 
tion in the Eo-ensemble we obtain the scattering function 
in the E-ensemble in homogeneous field as 



S-\q) = T + cq'^ - 2iT(3nal E^ 



(qET 

g2 



(18) 



The isotropic term on the right-hand side of is due 
to interactions of the induced dipoles with the fluctuat- 
ing electric field to order , while the anisotropic term is 
associated with the dipole-dipole interactions. The com- 
parison of the electric-field contribution with the term 
T = T — Tc va (fT5|) yields that the critical temperature 
in an electric field is shifted upwards, which is in accor- 
dance with the thermodynamic consideration for the van 
der Waals gas in the next paragraph. The density fluctu- 
ations become anisotropic in an electric field, such that 
the fluctuations with wave vectors transverse to the 
field strength will be enforced, while fluctuations with the 
wave vectors q|| parallel to E remain unchanged. The up- 
ward shift of Tc is due to the isotropic term, which is ab- 
sent in the Eg-ensemble. The anisotropic character of the 
density fluctuations described by Eq. has the con- 
sequence that the instability limit for fluctuations with 
wave vectors transverse to E will be reached at higher 
temperatures, which also means that an ordered state 
with the interface parallel to the electric fleld will be pre- 
ferred. Note that the anisotropic term in (|18p has the 
same form as the uniaxial anisotropy in magnetic sys- 
tems (23 |. and can change the upper critical dimension 
from dc = 4 to dc = 3. 

As a last point in this subsection we will check the 
validity of the relation between the compressibility and 
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the structure factor 

l + nJh{r)d^r^S{q^O)^kBT(^^^ (19) 

in an electric field in the E-ensemble at the critical point. 
The computation of {dp/dn)rp using Eq. (pS]) in the 
next Section yields at the critical point {dp/dn)rp — 
— ^alficE'^- The latter should be compared with 
S^^{q — > 0)\T=Ta which is obtained from Eq. (fT8|l as 

The average of the right-hand side over the directions of 
the unit vector = q'^/q gives the factor 6'^'^ /3, so that 
(fT9ll is fulfilled. 



obtains immediately from ((2T|) 



id 



,)<p?>+- 



(22) 



where Qp{E) — (47rsinh(x)/x) with x = PpE is the 
partition function of the non interacting dipole moments 
in an electric field, and < pi >= p (— 1/x + coth(x)) is 
the absolute value of the mean dipole moment. The cal- 
culation of occurs similar to that in ^ . The contri- 
bution to the free energy can be obtained from Eq. ^ 
by the replacement aEo — >< pi >=< p{E) >, and yields 



-At: N 
T{E) = F^+ Tj,{E) +N—< p{E) >^ - 



(23) 



B. Description of polar systems in the E-ensemble 

In this subsection we will generalize the statistical me- 
chanical description of polarizable liquids in an electric 
field developed above to polar liquid systems. As in- 
duced dipole moments are in many cases much smaller 
than the permanent ones, we neglect the polarization ef- 
fects. However, a generalization of the formalism taking 
into account both the effects of permanent and induced 
dipole moments is rather straightforward. In particular, 
the latter is expected to be of interest for the computa- 
tion of the dielectric constant for polar liquids. 

The interaction energy is now given by -ffez.d — Hdd — 
J^tPi^oiri) with 



,)P^ 



(20) 



being the energy of dipole-dipole interactions. The con- 
figurational integral in the E-ensemble is defined simi- 
larly to Eq. ©J where Hpoi in the expression of QAr(£'o) 
in Eq. ^ is replaced by H^i^d- The functional integration 
over Eo{r) results instead of in the following exact 
expression 



_ln(_detAo-i) 



QAr(E)= J drexp(-/3i/o)exp 
X / d^pi ■ ■■ d^pN 



X exp 



'/3Hdd+^ /dV(E(r)+47rP(r))' 
Svr J 



(21) 



where Ao{r — r') = S{r — r'), and P(r) = J2i Pi^i''" ~ ^i) 
is the density of the dipole moment. 

To enable a quantitative conclusion on the effect of 
dipole-dipole interactions we restrict ourselves here to 
the calculation of (3Ar(E) to the first order in Hdd- One 



wherejrp(£') — —kBT\nQp{E). The contribution of per- 
manent electric dipoles to the binary interaction energy 
in the leading order is given by 

W2,dd = IT.<p"> ^"^('■^ - <P^>- (24) 



The contribution of W2^dd to the Ornstein-Zernike for- 
mula can be obtained from ([Ml in a straightforward way, 
and is similar to Eq. mi). The last term in Eq. ((251) 
is similar to that for the free energy in the Eg-ensemble 
given by Eq. and results, as we will see in the next 
section, in a downward shift of the critical temperature. 



III. SHIFT OF THE CRITICAL 
TEMPERATURE IN AN ELECTRIC FIELD 

Using the above results we will consider the shift of the 
critical point of a simple liquid in an electric field. The 
pressure p — ~ {dT / dV)rp -^^ (not to be confused with 
the absolute value of the permanent dipole moment) is 
obtained from Eq. ^ as 



P = Po 



(25) 



where the pressure po in the reference state can be iden- 
tified with that of the van der Waals equation. The last 
term in (|25[) corresponds to the dipole-dipole interactions 
of induced dipoles, which are repulsive according to (|25|) . 
The interaction of the induced dipoles to the lowest or- 
der in as (= asEo) with an homogeneous field does not 
depend on volume, and thus does not contribute to the 
pressure. For a dielectric material in a parallel plate ca- 
pacitor this case corresponds to the charged capacitor in 
a open circuit. For computation of the pressure in the 

E-ensemble we use the relation p — ~ (j)T/dV^ with 
T ^ F + J'E'^{r)d^r/8Tr (the latter is the energy of the 
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external field), and F ~ -kBTlnQN^E) with Qn(E,) 
defined by Eq. ([5]), and obtain to order 



P^Po 



2tt 



2-rr,2 



(26) 



It follows from Eqs. (|25l26p that the shift of the crit- 
ical point in an electric field depends on the ensemble, 
i.e. if the electric circuit is open or closed in the exam- 
ple with parallel plate capacitor. The contribution of the 
electric field to p in Eq. ([25l) . which is due to the dipole- 
dipole interactions of induced dipoles, is positiv, and di- 
minishes the constant a of the van der Waals equation, 
and thus results in a decrease of the critical temperature 
Tc- It follows from Eq. ([26|l that the contribution of the 
electric field to the pressure at constant E increases the 
value of the constant a, so that the shift of the critical 
temperature is upward in the E-ensemble. Note that Tc 
for the van der Waals equation is given by the expres- 
sion fc^Tc = 8a/276. With the shift of a derived from 
Eq. ([261), Sue = {2n/3)al'E'^, we obtain the shift of the 
critical temperature as 6Tc = IGirpcO^E'^ /27kB- The 
latter coincides with the expression given by Landau and 

Lifshitz [3 (5T = ^ (0) /-^ if evaluated with the 

Clausius-Mossotti formula and the van der Waals equa- 
tion. The upward shift in the E-ensemble can be under- 
stood as follows. The macroscopic electric field is given 
by averaging Eq. ([2]) . The uniformity of E is ensured by 
fluctuations of the external field, which compensate the 
fluctuations of the electric held of induced dipoles 



SEo 



]T{r- 



which are due to the thermal motion of atoms. Thus, 
the interactions of the induced dipoles with the external 
field, ~ —piSEo with pi ~ aiEo and SEq given above, 
are negative, and are of the same order of magnitude 
as the dipole-dipole interactions. The first term on the 
right-hand side of the effective binary energy in Eq. (fT2|) 
originates from this term. It follows from Eq. (fT2|) or 
Eq. that the total sign of the effective binary interac- 
tions is determined by the first term. In the E-ensemble 
the interactions of the induced dipoles with the electric 
field enforce spontaneous inhomogeneities of the density, 
leading to an increase of the electric energy. This is sim- 
ilar to the behavior in an inhomogeneous electric field, 
which favors demixing 

The shift of the critical temperature temperature for 
polar liquids can be obtained from Eq. (|23p similar to 
the derivation of the Tc shift in the Eo-ensemble at the 
beginning of the Section IIIII The contribution to the 
pressure is obtained from Eq. ([25]) by the replacement 
aEo p{E) >. Thus, the permanent dipole moments 
result in a downward shift of the critical temperature. As 
the average dipole moment in ()23p equals for moderate 
field Ep^jS/i the downward shift of the critical tempera- 
ture is proportional to the square of the field strength. 



The results of this paper show that the sign of the shift 
of the critical temperature is determined by the perma- 
nent dipole moments of the atoms or molecules. Un- 
der the condition p/ aE > 1 the shift of Tc is downward 
for polar systems and upward for unpolar, but polariz- 
able systems. This conclusion is in accordance with the 
available experimental results in [3| , [IB] j [IB] j [HI , [IB] . 
The molecule sulfur hexafluoride {SFq) does not possess 
a permanent electric dipole moment, so that according to 
our work an upward shift is expected. Such an upward 
shift for SFq was established in [IB]. A downward shift 
for nitrobenzene {CqHzN02), which possesses a perma- 
nent dipole moment of 4.23 D, was measured in [l9|, and 
in the earlier work [3] , and is also in accordance with our 
prediction. A downward shift was also measured in [IB] 
for polystyrene/cyclohexane system. Styrene monomers 
have a dipole moment of 0.13 D, while the cyclohexane 
molecules are unpolar, so that the negative Tc shift in this 
system is also in accordance with our prediction. In the 
recent paper [IB] a downward shift of the disorder to or- 
der transition in polystyrene-block-polyisoprene diblock 
copolymers in concentrated solutions was measured. For 
flelds as high as 8.5 kV/mm a decrease in Tc by more 
than 1.5 K was observed. The dipole moments of styrene 
and isoprene monomers are equal to 0.13 D and 0.38 D, 
respectively, which explains the downward shift. 

The estimate of the electric energies of the isoprene 
monomers due to the polarization effect and the per- 
manent dipole moment p in IJ] yields for the ratio 
pI aE 3> 1, and thus legitimate the neglection of induced 
electric dipole moments in Section [II Bl 



IV. CONCLUSIONS 

To conclude, the statistical mechanical description of 
liquid systems in an electric field in the ensemble at con- 
stant potential, which we establish first for polarizable 
liquids, provides the basis for systematic calculations of 
thermodynamic quantities, dielectric constants, nonlin- 
ear effects in E beyond E^ for broad class of liquid sys- 
tems like mixtures, polymers, etc. The generalization of 
this formalism to liquid systems with permanent dipole 
moments enables us to analyze the shift of critical tem- 
perature in an electric field in polarizable and polar liq- 
uids. Our prediction that the shift is upward for po- 
larizable systems, and downward for polar liquids is in 
accordance with available experimental results. 



Acknowledgments 

Financial support from the Deutsche Forschungsge- 
meinschaft, SFB 418 is gratefully acknowledged. We 
would like to thank the Referees for bringing our atten- 
tion to t he p aper [25j . and A. Boker for sending us the 
preprint (25| prior to publication. 



6 



[1] J.G. Kirkwood, J. Chem. Phys. 4, 592 (1936). 

[2] R. Zwanzig, J. Chem. Phys. 25, 211 (1956). 

[3] A. Isihara and R.V. Hanks, J. Chem. Phys. 36, 433 
(1962); A. Isihara, J. Chem. Phys. 38, 2437 (1963). 

[4] L.D. Landau and E.M. Lifshitz, Electrodynamics of Con- 
tinuous Media, Pergamon Press, Oxford, 1984. §18. 

[5] Y. Tsori, F. Tournilhac, and L. Leibler, Nature 430, 544 
(2004). 

[6] K. Amundson et al., Macromolecules 26, 2698 (1993). 
[7] T.P. Russell et al. , Macromolecules 23, 890 (1990). 
[8] E. Gurovich, Phys. Rev. Lett., 74, 482 (1995). 
[9] T. Thurn-Albrecht, J. DcRouchey, T. P. Russell, Macro- 
molecules 33, 3250 (2000). 
[10] Y. Tsori, D. Andelman, C.-Y. Lin, M. Schick, Macro- 
molecules 39, 289 (2006). 
[11] Y. Tsori and D. Andelman, J. Polym. Scie. Part B - 

Polym. Phys. 44, 2725 (2006). 
[12] K. Schmidt et al. , H.G. Schoberth, M. Ruppel, H. Zettl, 
H. Hansel, T.M. Weiss, V. Urban, G. Krausch, and A. 
Boker, Nature Materials 7, 142 (2008). 
[13] M.W. Matsen, J. Chem. Phys. 124, 074906-1 (2006). 
[14] P. Debye and K. Kleboth, J. Chem. Phys. 42, 3155 



(1965). 

[15] D. Beaglcholc, J. Chem. Phys. 74, 5251 (1981). 

[16] D. Wirtz and G.G. Fuller, Phys. Rev. Lett. 71, 2236 
(1993); D. Wirtz, K. Berend, and G.G. Puller, Macro- 
molecules 25, 7234 (1992). 

[17] A. Ormki, Europhys. Lett. 29, 611 (1995). 

[18] G.A. Zimmerli et al., Phys. Rev. Lett. 82, 5253 (1999). 

[19] K. Orzechowski, Chem. Phys. 240, 275 (1999). 

[20] J. Hegseth and K. Amara, Phys. Rev. Lett. 93, 057402-1 
(2004). 

[21] C. Park et al.. Polymer 44, 6725 (2003). 

[22] R. Kubo, Statistical Mechanics, North-Holland Publish- 
ing Company, Amsterdam, 1965. 

[23] G.S. Rushbrooke, in Physics of Simple Liquids, ed. by 
H.N.V. Temperley, J.S. Rowlinson, and G.S. Rushbrooke, 
North-Holland Publishing Company, Amsterdam, 1968. 

[24] A. Aharony, in Phase Transitions in Critical Phenomena, 
Vol. 6 (Academic Press) 1976. 

[25] H.G. Schoberth, K. Schmidt, K. Schindler, and A. Boker, 
unpublished. 



